Ratchet potential and rectification effect in Majorana fermion SQUID 
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Motivated by a recent experimental progress in realizing Majorana fermions (MFs) in 
a heterostructure of a spin-orbit coupling nanowire and superconductor (V. Mourik et al. 
Science. 1222360), we investigate a SQUID formed by the novel superconductor-nanowire- 
superconductor Josephson junction which contains MFs and a conventional superconductor- 
insulator-superconductor junction. It is shown that the critical current of the SQUID is different for 
the two current directions. Since the asymmetric Josepshon current forms a ratchet potential for 
the dynamics of superconducting phase, a rectification effect is expected when the SQUID is driven 
by an ac current. These novel properties are expected to be useful for probing the elusive MFs as 
well as for their dynamics. 

PACS numbers: 74.50.+r, 85.25.Dq, 71.10.Pm 



Introduction. - Majorana fermions (MFs) in topological 
superconductors have drawn much attention recently^— , 
due to their potential application in topological quantum 
computatio n 1 1 i 16 ~ —. They have been predicted to ex- 
ist in a number of systems, including the spin-triplet p- 
wave superconductor— superconductor-topological in- 
sulator interface-, and semiconductor-superconductor 
heterostructur o 7 ' 9 ' 10 . 

Among all these candidate systems, the one dimen- 
sional nanowire-superconductor heterostructure is of par- 
ticular interest* 9 -*^, due to the well developed nanowire 
manufacturing technique. It is composed of a spin-orbit 
coupling nanowire in proximity to an s-wave supercon- 
ductor under a moderate magnetic field, and MFs are 
expected at the ends of the wirei. Actually, in a re- 
cent experiment, a nanowire-superconductor device was 
fabricated with an InSb nanowire connecting to super- 
conductor NbTiNSS. A zero-bias peak in the differential 
conductance has been detected, in agreement with the 
theoretical prediction of the resonant Andreev effect^, 
which is considered as a strong evidence for the existence 
of MFs. 

Here we propose to build a SQUID using the 
MF nanowire Josephson junction successfully fabri- 
cated in this recent experiment^ and a conventional 
superconductor-insulator-superconductor (SIS) Joseph- 
son junction, as schematically depicted in Fig. 1. We 
predict a directionally asymmetric critical current in this 
SQUID due to the unconventional current phase relation 
(CPR) involving MFs. This direction-dependent criti- 
cal current is ubiquitous for both adiabatic and fast MF 
dynamics thus can provide a supportive evidence for the 
MF existence. Since this asymmetric critical current cor- 
responds to a ratchet potential, i.e. a periodic potential 
without reflection symmetry, for the "particle" of super- 
conducting phased, there appears a novel rectification 
effect in the SQUID, that is, an ac driving current in- 
duces a rectified dc voltage. This rectification effect is 
not only dominated by the dynamics of phase differences 
in Josephson junctions, but also strongly influenced by 




FIG. 1: Schematic setup of SQUID composed of a novel junc- 
tion containing two Majorana fermions and a conventional 
SIS junction. 



the dynamics of the MF state. 

Current-phase relation. - The MF Josephson junction 
reported in the recent experiment consists of an InSb 
nanowire in between two NbTiN superconductors, with 
a magnetic field along the nanowire; two MFs are ex- 
pected at the two sides of the tunneling barrier induced 
by the gate voltage at the middle of the nanowire^. This 
MF Josephson junction can be described by a tunneling 
Hamiltonian H = Hq+Ht- Here Ho is the Kitaev model 
for two semi- infinite topological superconducting wires 1 , 

H =H L + H R (1) 
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where |A| is the superconducting gap, 7^ = e~ l ~Cj + 
e l ~Cj, jj t B = i(e~ l ~ clj — e l ~Cj) are the Majorana op- 
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erators, with c] the spinless electron creation operator, 
tfj = <pl for j < 0, and = for j > 1 are the super- 
conducting phases of the two superconductors. The two 
end MFs 70,5 and 71 ,a do not appear in the above Hamil- 
tonian, which implies a two-fold degenerate ground state 
for Hq. The two superconductors are connected through 
the electron tunneling, 



Ht = T e ct cr + h.c. 



(2) 



where cj, = cq, cr — c\, with T e the electron tunneling 
matrix. 

The Josephson current can be evaluated by the charge 
change in the left superconductor, 
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where T — T e /2, Nl is the electron number operator of 
the left superconductor, 9 = <p R — <pl is the phase dif- 
ference between the right and left superconductors, and 
\G) is the wave function of the ground state. 

For T << |A|, the Josephson current can be calcu- 
lated by the perturbation method. Since the ground state 
of the unperturbed Hamiltonian Hq is two-fold degener- 
ate, to the lowest-order perturbation, the wave function 
\G) = I Go) is contributed merely from the ground state 
subspace of Hq. Therefore, Eq. (3) reads, 



I = J M sm(9/2)(Go\i~fLjR\G ) 



(4) 



with 7l = 7o,B) Jr = 7i,Aj an d the critical current Jm — 
eT/h. The lowest-order perturbed state |Gq) evolves with 
time according to the Schrodinger equation, 

ihj t \G ) = {Tco S (6/2)(-i lLlR )} \G Q ), (5) 

where only one term in Ht expressed by MF operators 
has been registered, since other three terms project \Gq) 
out of the subspace of degenerate ground states of Hq. 

In more general cases, the end MFs may slightly devi- 
ate from the zero energy due to the correlation to other 
parts of the superconductors they attached&i^. This 
small coupling can be taken into account generally by 
introducing two terms 8l"/l + 8r~/r into the Hamilto- 
nian Hq, with 8l and S R including all possible interac- 
tions to 7l and j R respectively. As far as MF behav- 
iors are observed, they should be small, and we take 
5 1, = 5r = 6 « T in the present work. Since these 
two terms account for possible interactions inside each of 
the two superconductors, they do not change the current 
expression (4). The dynamic of the end MFs is however 
enriched, and Eq. (5) is modified into, 



iHj t \Go) = [Tcos(0/2)H7L7fl) + KlL + 7fl)] \G ). (6) 



The ground state of H Q can be described by the two 
eigenstates of the MF operators —ijLj R \±.) — ±|±) as 
\G Q ) = Pi|+) + P 2 \-}. In this basis, Eq. (6) takes a 2x2 
matrix form, 



n d ( P i 
lH dt P2 



Tcos(6»/2) S-iS 
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Pi 
Pi 



,(7) 



and the Josephson current is 

lM(d,Pi,P 2 ) - Jm(|P 2 | 2 - |Pi| 2 )sin(0/2). (8) 

We notice that if the coupling S equals zero, Eq. (7) is 
reduced back to the diagonal one in Eq. (5), and eigen- 
states (Pi,P2) are (1,0) or (0,1) corresponding to the 
two different parities. This leads to a 47r-period CPR 
known as the fractional Josephson effect^ ' 14 ! 23 , 



I M (9) =±J M sin(0/2), 



(9) 



as seen in Fig. 2a. With the coupling term 8, the par- 
ity is no longer a conserved quantity, and the evolution 
of I Go) becomes a Landau-Zener problem 2 ^, which has 
been studied in many other systems and exhibits rich 
dynamics. There are two limits in the dynamics gov- 
erned by Eq. (7): For the fast process in which the phase 
difference changes quickly, the coupling 8 has no much 
influence and the system behaves essentially in the same 
way as for 8 — 0. However, for the slow adiabatic pro- 
cess, the instant ground state of the 2x2 matrix is always 
reached, and the Josephson current is given by diagonal- 
izing Eq. (7), 



I M {0) = J M A(0)8m(0/2), 



(10) 



where 
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As seen in Fig. 2a, the CPR is of 2ir period in this 
caseS&ii. It is clear from Eq. (8) that the CPR of the MF 
Josephson junction strongly depends on the dynamics of 
the end MFs, and generally deviates from the simple si- 
nusoidal function of phase difference realized in conven- 
tional Josephson junctions. 

Asymmetric critical current. - We connect the MF 
Josephson junction to a SIS Josephson junction, and form 
a SQUID schematically displayed in Fig. 1. The CPR of 
the SIS junction is the conventional one In — Jn sin 9n, 
where Jm is the critical current and 9n is the phase dif- 
ference. With a magnetic flux $ in the SQUID, the 
phase differences in the two junctions are related by 
9n = 9 + 27r$/<I>o, where $ = h/2e is the flux quan- 
tum. The Josephson current of the SQUID is 



1(9) = Jjq sin 



2tt$ 



Im(0). 



(12) 
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FIG. 2: (Color online), (a) CPR of MF Josephson junc- 
tion for the fast process (solid line) and adiabatic process 
(dotted line). CPR of the SQUID for the fast process (b) 
and adiabatic process (c) with typical applied magnetic flux 
$ = (dotted line), $o/8 (dashed line), and $o/4 (solid line), 
(d) Interference pattern for the critical current of the SQUID 
in the fast process (solid line) and adiabatic process (dotted 
line). Parameters are taken as 5 = 0.01T and Jjv = Jm- 



We plug Eqs. (9) and (10) into the above equation, and 
obtain the CPR of the SQUID for the fast process and 
adiabatic process respectively. The results for three typ- 
ical $ values are shown in Figs. 2b and 2c. 

It is interesting to observe that, for a finite magnetic 
flux, the currents in the opposite directions, i.e. I > 
and I < 0, arc different in both fast and adiabatic pro- 
cesses. As an experimentally observable quantity, we 
evaluate the magnetic flux dependence of the critical 
Josephson current. As displayed in Fig. 2d, the critical 
Josephson currents are different in the two opposite di- 
rections for a given magnetic flux, in contrary to the well- 
known interference pattern in a conventional SQUID. 

The asymmetric CPR is caused by the sub-harmonic 
function sin(0/2) in Eq. (9) for the fast process, and the 
non-sine function in Eq. (10) for the adiabatic process, 
both originated from the existence of MFs as discussed 
above. It is due to point out that a SQUID formed by 
two conventional Josephson junctions cannot realize this 
asymmetric CPR even if they possess different critical 
currents. Observed both for fast and adiabatic dynamics, 
the asymmetry of critical Josephson current with respect 
to the flowing direction is ubiquitous for MF systems, as 
compared with the fractional Josephson effect which is 



expected only for fast dynamics. 

The critical current asymmetry in this SQUID struc- 
ture can serve as an ideal probe for the end MFs. From 
the experimental perspective, the MF Josephson junc- 
tion has already been fabricated and its I — V curve 
measured 2 ^. Our setup can be built simply by connecting 
it to an additional SIS junction and forming a SQUID, 
which is not a difficult task in experiments. Meanwhile, 
the measurement of critical current in junction systems 
is already a routine procedure, leaving no difficulty in 
experimentally detecting the directional asymmetry. 

Rectification effect. - The asymmetry in critical cur- 
rents for two flowing directions in the MF SQUID re- 
vealed above has an interesting consequence known as 
the rectification effect under an ac driving current, which 
has been discussed in various contexts so fa r 22 i 25 ~ -2£. In 
the present MF system, the "particle" of superconduct- 
ing phase feels a ratchet potential as seen in Fig. 2 due to 
the unconventional CPRs in Eqs. (9) and/or (10), where 
the translational anti-symmetry 1(9 + tt) = —1(6) en- 
joyed by conventional SQUIDs is broken 28 . When the 
amplitude of the ac driving current is chosen in between 
the two different critical currents, pure supercurrent flow 
is possible only in one direction, while the current in the 
opposite direction exceeds the critical current and con- 
tains a normal part. In this rectified voltage is 
induced across the junction. As the degree of asymmetry 
in critical current is different for the adiabatic and fast 
processes, this rectification effect is strongly influenced 
by the MF dynamics. 

To be specific, we adopt the resistively and capaci- 
tively shunted junction (RCSJ) model, with the phase 
difference 9 evolves with time according to the dynamic 
equation, 



hC d 2 9 
1 2e~~dl I 



h d9 



2eR dt 



-I(t) 



- Jn sin 



2tt$ 



(13) 

J M (\P2\ 2 -\Pi\ 2 )sm(9/2), 



where I(t) is the external current, and R and C are the 
resistance and capacitance of the SQUID. The dynam- 
ics of the SQUID can be obtained by solving Eqs. (7) 
and (13) simultaneously. The resistance of the SQUID 
is taken as R = 2kf2, in the same order of that for the 
MF Josephson junction in the recent experiment 2 ^. As 
a small hysteresis loop in the I — V characteristics was 
observed experimentally 2 ^, we focus here on a slightly 
underdamped case with Q-factor Q = eRy / 2CT/h = y/2. 

Without loosing generality, let us put Jn = Jm, and 
take $ = $0/3 as an example, where for the adiabatic 
process the left and right flowing critical currents are 
0.866Jm and — 1.77Jm respectively, while for the fast 
process, they are 1. 37 Jm and — 1.97Jm respectively, as 
seen in Fig. 2d. 

Now we apply an ac driving current I(t) = Je sin(uit), 
with Je = Jm- The time evolutions of the phase dif- 
ference 9(t) and the MF state \Go(t)) are shown in Fig. 
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FIG. 3: (Color online). Dynamics of the phase difference 9(t) 
(a,b), and the MF state \Go(t)) (c,d), under an ac driving 
current Je = Jm, with driven current frequency ftw = O.Ol7r<5 
(a,c), and fiu = Q.\ti5 (b,d). Other parameters are taken as 
Jjv = Jm, R = 2fcfi, Q = V2 ( see text for detailed discus- 
sions), and 8 = 0.01T. 
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FIG. 4: (Color online). Same as in Fig. 3 except for Je ~ 
1.5 Jm • 



3, for two typical frequencies hu = 0.0±7r<5 and 0.l7r<5. 
It is clear that for hw = O.Ol7r<5, the MF state evolves 
away from the initial state and takes the instant ground 



state, a typical adiabatic behavior. Simultaneously, we 
find that the phase difference increases with time mono- 
tonically, which indicates a dc voltage. This is the voltage 
rectification effect expected from the asymmetric critical 
current. For haj = 0.l7T<5, however, the MF state stays at 
the initial state, a typical fast non-adiabatic process be- 
havior, and the phase difference is oscillating around zero 
indicating no rectified voltage in the system. The differ- 
ent rectification behaviors for different driving current 
frequencies come from the value Je , which is in between 
the two critical currents of different directions for the 
adiabatic limit, but below the two critical current for the 
fast process. For Je = 1.5Jm, which is in between the 
two critical currents in both limits, the time evolution of 
the phase difference and the MF state is shown in Fig. 4. 
We observe the rectification effect for both frequencies. 
With Jm — 4:0nA suggested from the experiment^, the 
rectified voltage V — (8)h/2e = 0.2HJm /2e 2 is in order 
oflO/iV. 

We numerically simulate various junction parameters, 
and find this rectification in a large parameter space. The 
experimental observation of this rectification effect is ob- 
viously helpful for identify the existence of MFs. Mean- 
while, since it disappears when the driven frequency ap- 
proaches S under appropriate driving currents, it is pos- 
sible to draw information on the MF coupling. The rec- 
tification effect is scalable, and the rectified voltage can 
be amplified simply by connecting several SQUIDs in se- 
ries, making its experimental observation and potential 
application plausible. 

Conclusion. - In summary, we have demonstrated 
that Majorana fcrmions produce unconventional current- 
phase relations in a Josephon junction including a spin- 
orbital coupling nanowire, with the detailed forms de- 
pending on the dynamics of the Majorana fermions. The 
SQUID formed by this novel Josephson junction and a 
conventional one shows ubiquitously direction-dependent 
critical currents. This ratchet potential established by 
Majorana fermions for the dynamics of superconducting 
phase yields a rectified dc voltage across the SQUID when 
it is driven by an ac current. These novel properties 
can be explored as useful probes to the elusive Majorana 
fermions. 
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